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steady aerodynamic influence coefficient matrix (eq. 3.5-41, 
ref. 1) 

= [T^] [q c ] [A pQ ] 1 (eqs . 5.3-3, ref. 1) 
semispan of a thin body 
pressure coefficient 

matrix of isolated thickness induced pressure coefficient 

matrix of interference coefficients relevant to vortex 
interference 
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matrix of leading-edge panel force per unit panel span 

components of aerodynamic force along Body Axes 

matrix of aerodynamic force components induced by thickness 
but not including thickness interference forces 

matrix of interference aerodynamic force components induced 
by thickness 

components of moment of inertia about Body Axes 
identity matrix 
Mach number 

rolling, pitching, and yawing moments about Body Axes 

matrix of Mach number variation for calculating the effect of 
isolated thickness on body forces; 
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P 
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rate of roll 

= ~ , nondimensionalized P 
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structural transformation matrix — from nodal displacement 
components to aerodynamic control point rotations 
(eq. 4.2-101, ref. 1) 

Q 

rate of pitch 

Q 

~ 2U 9 nondimensionalized Q 

q 

dynamic pressure apparent to Reference Axis System 
(eqs. 2.3-75, ref. 1) 

^00 

dynamic pressure apparent to the Fluid Axis System 
(eqs. 2.3-74, ref. 1) 

bq t 4 

dynamic pressure distribution formed as [ q ] but does not 
include the coordinates of the slender body mean inter- 
ference surface panels 


dynamic pressure distribution referring to the aerodynamic 
centroids of all aerodynamic segments (i.e., slender body 
centerline segments, slender body mean interference sur- 
face panels, and thin body mean surface panels) 

(eq. 3.5-54, ref. 1) 

R 

rate of yaw 

R 

= -Tpj , nondimensionalized R 

S w 

wing area 

[SD S ] 

transformation matrix — from a Stability Axis System to a 
Body Axis System (eqs. 5.3-44, ref. 1) 

U 

X-component of the translational velocity of the Body Axis 
System relative to the Inertial Axis System 

X,Y,Z 

Reference Axis System (fig. 1) 

^ * Y B ’ Z B 

Body Axis System (fig. 1) 

NY c s! 

matrix of the coordinates of the aerodynamic centroids of all 
aerodynamic segments (i.e., slender body centerline seg- 
ments, slender body mean interference surface panels, and 
thin body mean surface) 

b? t 4 

matrix of the coordinates of the aerodynamic centroids of all 
aerodynamic segments excluding the coordinates of the 
slender body mean interference surface panels 
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angle of attack 
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Subscripts : 

B 

E 

E* 

I 
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sideslip angle 
flight path angle 
eleven angle 
aileron angle 

elevator angle (6 is substituted in place of 6 e in the 
longitudinal stability computation) 

rudder angle 

pitch attitude 

matrix of elastic rotations at control points 
angle of bank 

matrix of flow incidence due to angle of attack 

matrix of flow incidence due to sideslip 

matrix of flow incidence due to aileron deflection 

matrix of flow incidence due to elevator deflection 

matrix of flow incidence due to rudder deflection 

matrix of flow incidence arising from camber shape 

matrix of flow incidence arising from aircraft motion and 

control surface deflections 

matrix of flow incidence due to roll rate 

matrix of flow incidence due to pitch rate 

matrix of flow incidence due to yaw rate 
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quantity relevant to a flexible aircraft 
quantity associated with elastic increment 
Ith thin body 
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°° quantity referred to free stream 
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STATIC AND DYNAMIC STABILITY ANALYSIS OF 


THE SPACE SHUTTLE VEHICLE-ORBITER 
Wei J. Chyu, Ralph K. Cavin,* and Larry L. Erickson 
Ames Research Center 


SUMMARY 


The longitudinal static and dynamic stability of a Space Shuttle Vehicle- 
Orbiter (SSV Orbiter) model is analyzed using the FLEXSTAB computer program. 
Nonlinear effects are accounted for by application of a correction technique 
in the FLEXSTAB system; the technique incorporates experimental force and 
pressure data into the linear aerodynamic theory. A flexible Orbiter model is 
treated in the static stability analysis for the flight conditions of Mach 
number 0.9 for rectilinear flight (1 g) and for a pull-up maneuver (2.5 g) at 
an altitude of 15.24 km. Static stability parameters and structural deforma- 
tions of the Orbiter are calculated at trim conditions. For the dynamic sta- 
bility analysis, the characteristics of damping in pitch are investigated for 
a Mach number range of 0.3 to 1.2. The calculated results for both the static 
and dynamic stabilities are compared with the available experimental data. 
These comparisons show that it is necessary to use the correction technique 
due to the nonlinearity. The calculated results indicate that the elevon 
effectiveness and the longitudinal stability are lower for the elastic Orbiter 
model than for the rigid model. The g effect on the trim solutions is signif- 
icant for both the rigid and flexible models. Calculated pitch damping is in 
good agreement with the experimental results. 


INTRODUCTION 


Aeroelastic effects were neglected in most aircraft stability analyses 
until recently when the large transonic and supersonic aircraft posed diffi- 
cult stability problems. The apparent static stability characteristics of 
these aircraft are drastically affected by the elastic deformations of their 
low-aspect ratio thin wings and slender fuselages. Further, the dynamic sta- 
bility is significantly affected by the dynamics of the structure through 
unsteady aerodynamic coupling between the rigid body and the structural 
motions. For these aircraft, a mathematical model based fundamentally on the 
dynamics of a flexible body is required to predict static and dynamic stabil- 
ity. To meet this demand, a system of digital computer programs known as the 
FLEXSTAB system was recently developed (ref. 1). 

The distinctive features of the FLEXSTAB system are its ability to eval- 
uate static and dynamic stability, trim state, inertial and aerodynamic 
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loading, and elastic deformations of aircraft configurations at subsonic and 
supersonic speeds. The FLEXSTAB system was used previously in the analysis of 
such modern aircraft configurations as the Boeing 707 (ref. 2) and the YF-12A 
(ref. 3). Because the FLEXSTAB system imposes a restriction of small dis- 
turbances to airflow by the configurations, the configuration of the Boeing 
707 airplane, which falls within this restriction, results in good agreement 
between the theory and the experiments for the pressure distribution on the 
wing. However, a complex configuration like that of the YF-12A aircraft 
results in significant discrepancies between the theoretical results, wind- 
tunnel data, and flight- test data, particularly in the region of the forebody 
chine and the outboard wing. The discrepancies are attributed to the vortex 
flow in the forebody chine, a flow that the linear aerodynamic theory within 
the FLEXSTAB cannot simulate. 

The purpose of this report is to describe the technique of applying the 
FLEXSTAB system to the static and dynamic stability analyses of a flexible 
Space Shuttle Orbiter configuration. The SSV Orbiter configurations feature a 
bulky fuselage blended with a swept thick wing and forebody chine (ref. 4) . 
This type of configuration, as that of the YF-12A, is likely to cause flow 
separation on most of its surface and vortex flows emanating from the leading 
edge of the wing. The exact linear theory alone cannot treat this kind of 
nonlinear flow problem. The correction approaches available within FLEXSTAB 
are thus used in the static stability analysis by incorporating the experi- 
mental aerodynamic data into the linear aerodynamic theory. The linear aero- 
dynamic and structural theories are used with the wind-tunnel results for the 
rigid model to predict the aeroelastic effects of the actual airframe config- 
uration, and thus to give improved predictions of the actual aircraft stabil- 
ity and elastic deformation. For the dynamic analysis, the characteristics of 
damping in pitch are investigated using the linear theory of FLEXSTAB. 

For the static stability analysis the following analytical methods in the 
FLEXSTAB are used: 

1. Linear aerodynamic and structural theories — linear theory 

2. Implementation of experimental aerodynamic forces and moments into 
the linear- theories in (1) but using program-computed pressures — ■ force 
correction 

3. Implementation of experimental pressures into the linear- theories in 
(1) but using program- computed forces and moments (based on the experimental 
pressures) — pressure correction 

4. Implementation of experimental aerodynamic forces, moments, and lift- 
ing pressures into the linear theories in (1) — force and pressure correction 

The simple approach is to use all linear theory as in Method 1. When 
this was done it was discovered that the FLEXSTAB model did not always give 
good results for the pressure distributions or for the force and moment coef- 
ficients, the zero alpha coefficients being particularly poor. In addition, 
at the relatively high angle of attack at which trim occurs, the wind-tunnel 


2 



data exhibit a nonlinear variation with angle of attack. As mentioned above 
this nonlinear behavior cannot be predicted by the linear theory of FLEXSTAB. 
Consequently, Methods 2-4 were tried in an attempt to improve the linear 
FLEXSTAB results. 

As shown in reference 1 (vol. I, pp . 5-21) and also in the appendix, the 
flexible aircraft aerodynamic derivatives in the linear theory analysis are 
composed of the sum of two quantities — the rigid aircraft aerodynamic deriva- 
tives and an increment due to static aeroelasticity . In Methods 2-4, which 
utilize nonlinear wind-tunnel force and pressures, this separation into a 
rigid component and an aero-elastic increment is still assumed to be valid. 

Because of the discrepancies between the nonlinear experimental data and 
the linear theoretical results of Method 1, this method cannot be expected to 
predict an accurate trim solution for either the rigid or elastic case. Since 
the trim solution for the rigid case depends only on the integrated effects of 
the pressures (i.e., the pressure distribution can be wrong so long as the 
forces and moments are correct). Method 2, which replaces the theoretically 
computed forces and moments with nonlinear wind-tunnel values, is appropriate 
for the trim solution of a rigidly modeled vehicle. Since the aeroelastic 
increments in the flexible vehicle analysis depend explicitly on the pressure 
distribution. Method 2 would not be appropriate for the flexible aircraft 
analysis unless the pressure distribution is correctly predicted by the 
method. The incorrect pressure distribution that results from using 
Methods 1-2 thus will in general produce incorrect aeroelastic increments, 
even though correct wind-tunnel values are used for the rigid components of 
the force and moment coefficients. Thus, Method 3 was tried so that the cal- 
culated aeroelastic increments would be based on the measured pressures of the 
rigid wind-tunnel model rather than on the incorrect theoretical pressures by 
Methods 1-2. 

It was expected that the rigid components of the aerodynamic force and 
moment coefficients computed by FLEXSTAB from these experimental pressures 
would be in good agreement with the wind-tunnel results. Such was not the 
case, however, and there is no satisfactory explanation for the disagreement. 
Consequently, Method 4 was used. In Method 4, the experimental pressures are 
used for computing the aeroelastic increments, and the wind-tunnel force and 
moment coefficients were used (in place of the rigid component values 
FLEXSTAB computes from the experimental pressures) for computing the rigid 
aerodynamic derivatives and trim solution. 

It should be noted that the effect of the pressure correction in Method 3 
is to make corrections to the theoretical pressure terms which are present in 
both the rigid and elastic components of the theoretical aerodynamic force. 

As both the force and pressure are introduced in the correction (Method 4) , 
the effect of pressure correction is only on the elastic increment as its 
effect on the rigid component is superseded by the force correction. 

The main context of this report presents the techniques of applying the 
FLEXSTAB to the Space Shuttle Orbiter configuration. The analytical descrip- 
tion of the technique is summarized in the appendix which is based on the 
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voluminous FLEXSTAB program documents. The correction method is also derived 
in the appendix since there is a lack of clarity in the documents on the 
correction method. 

The computations for the longitudinal static stability study are made for 
two flight conditions: a' rectilinear flight (1 g) and a pull-up maneuver 

(2.5 g) , both at Mach number 0.9 and an altitude of 15.24 km. This flight 
condition is specifically chosen for the static analysis, because at this con- 
dition the greatest aeroelastic effects are observed within the Mach numbers 
investigated. In the dynamic stability analysis, a rigid Orbiter is treated 
using linear aerodynamic theory (ref. 1). The characteristics of damping in 
pitch, as given by the parameters and are investigated for a Mach 

number range of 0.3 to 1.2. Throughout this report, the analytical results 
are compared with available experimental data. 


METHOD OF ANALYSIS 


Geometry Definition of the Orbiter 

The Orbiter configuration (ref. 4) is described by using three types of 
component bodies in the Geometry Definition Program of the FLEXSTAB: (1) a 

slender body, (2) a thin body, and (3) an interference body (fig. 1). The 
Orbiter fuselage is represented by a slender body of revolution whose cross- 
sectional areas are the same as those of the fuselage. The line singularities 
representing the fuselage are extended downstream from the physical fuselage 
(order of 100 fuselage lengths) . This procedure was adopted to more accu- 
rately represent the actual flow field aft of the blunt-base fuselage (ref. 1, 
pp. 3-101, vol. IV, and/or vol. II). The camber line for the slender body is 
the locus of the area centroids of fuselage. The wing is represented by a 
thin body which is essentially the projection of the wing planform onto a mean 
plane. This projected wing plane is subdivided into panels with which 
FLEXSTAB associates appropriate aerodynamic singularities and control points 
(fig. 1). The 3.5° dihedral of the wing was neglected in order to ensure com- 
patibility with structural data. Wing camber and thickness effects are 
described by the wing airfoils (fig. 2) at various spanwise stations. The 
blunt airfoil trailing edge of the Orbiter has been replaced with a sharp 
trailing edge in keeping with a FLEXSTAB requirement. The Orbiter vertical 
tail is excluded in the geometry description. Studies of the rigid case indi- 
cate that the omission of the tail and the wind dihedral angle in the geometry 
definition has little effect on the results of the longitudinal stability 
analysis. Interference effects due to wing- induced downwash on the fuselage 
are approximately accounted for by introducing an interference body of polyg- 
onal cross section (fig. 1). The vortex singularities associated with the 
panels on the interference shell are used to account for the normal component 
of flow induced by the wing. In defining the Orbiter geometry the body flap 
was not incorporated into the present analysis because the flap would inter- 
sect the interference body in a manner which cannot be accommodated by the 
existing FLEXSTAB. 


4 



A rather lengthy sequence of computer runs was conducted for the rigid 
vehicle in an effort to determine what combinations of geometric parameters 
result in aerodynamic stability data that more closely correlate with avail- 
able wind-tunnel data* The effort involves changing shape and locations of 
the interference body, and modifying the paneling scheme and singularity dis- 
tribution to obtain an improved estimate for lift and moment-curve slopes. It 
results in the Orbiter geometry shown in figure 1, which consists of 50 con- 
trol points equally spaced along the fuselage, 98 finite-element panels on the 
wing, and 120 panels on the interference body* Computations with a larger 
number of finite-element panels and more densely defined airfoil than that 
shown in figures 1 and 2 were also tried but they yielded no further conver- 
gence in the numerical results . 


The Structural Model 

The External Structural Influence Coefficient (ESIC) program was used to 
couple FLEXSTAB to externally generated structural and inertia representations 
for the Orbiter. The structural model was provided by the Structures Group of 
Rockwell International Space Division. The model was a reduction of a large 
finite-element model of the Orbiter existing in Rockwell’s computer program, 
ASKA. An interface was created between the output data tape obtained from 
ASKA and subprogram ESIC of the FLEXSTAB system. 1 

The model provided by Rockwell International has a very dense nodal 
representation for the wing; it contains 200 nodes, each of which possesses 
the Z degree of freedom (fig. 3). In the leading-edge area of the wing, where 
no structural nodes are available, the FLEXSTAB interpolates the structural 
properties of nearby nodes to the aerodynamic centroid of the proper aero- 
dynamic model. 

The representation of the fuselage structure involves difficulties in 
representing, by a single line of nodes located along the centerline of the 
fuselage, the structural asymmetry induced in the fuselage by the upper "clam 
shell" cargo doors. The fuselage structural model provided by Rockwell Inter- 
national contains nine structural nodes each possessing X and Z degrees of 
freedom. The necessary structural data for the rigid Orbiter computation are 
shown in figure 1, and are input directly to the FLEXSTAB program. 


Aerodynamic Representation 

Static stability analysis - One of the major concerns during the course of 
this study was to select an aerodynamic representation which, together with an 
appropriate Orbiter geometry definition, accounts for nonlinear flow effects 
on the aerodynamic stability. It is evident that the nonlinear flow phenom- 
ena, such as the flow separation on most of the Orbiter surface and the vortex 
flows on the swept wing, cannot be described mathematically by the linear 
aerodynamic theory within the FLEXSTAB. However, the linear theory solutions 


^he interface was prepared by A. R. Curtis, Lockheed-California Company. 
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are of value in providing general aerodynamic trends; they also provide 
theoretical reference to which a correction on the deficiency of the theory 
can be made. A sequence of computations was conducted applying the following 
methods of aerodynamic representation: 

1. Linear aerodynamic and structural theories — linear theory 

2. Implementation of experimental aerodynamic forces and moments into 
the linear theory solution of (1) but using program computed pressures — 
force correction 

3. Implementation of experimental lifting pressures into the linear- 
theories solution of (1) but using program- computed forces and moments (based 
on the experimental pressures) — pressure correction 

4. Implementation of experimental aerodynamic forces, moments and lift- 
ing pressures into the linear theory solution of (1) — force and pressure 
correction 


The analytical details of the above methods are described in the appendix 
of this report; thus only the summary of the correction techniques and the 
computational procedures are described in this section. The experimental 
forces and moments (Cl, Cq, C m ) data are put into tabular forms as a function 
of angle of attack and elevon deflection, as shown in figure 4, and are pro- 
vided as inputs to the FLEXSTAB. The experimental lifting pressures are 
implemented into the FLEXSTAB by representing the pressures on both the wing 
and fuselage by a linear function of the angle of attack; that is. 


WT WT 

AC p. = AC p. 
i i 


+ AC 


WT 


( 1 ) 


where the subscript i denotes the ith body station, and superscript WT 

denotes wind-tunnel data. The coefficients AC^J T and AC^ are chosen so 

F i F i 

o a 

that the nonlinear experimental data for 0° < a < 14.5° are approximated by 
the analytically fitted expression (1) in the range of selected parameters; 
namely, a in the present consideration. The choice of this range of a is 
justified if the calculated trim angles fall within the range. The optimal 
linear approximations, as those shown in -figure 5, are applied to 50 fuselage 
stations and to 98 wing stations. 

In the approach for implementing the experimental data into the linear 
theory, the separation of the aerodynamic derivatives into a rigid component 
and into an aeroelastic increment, as in the linear theory analysis, is still 
assumed to be valid; that is. 
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( 2 ) 


fF 

^A 


V 


AF^ 1 

“a 

F 

zb a 

= 

Fzb j 

+ 
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zb a 



X 


am yb 
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where A. represents a, 6 e , or Q in the case of longitudinal analysis. The 
effect of the experimental aerodynamic force and moment implementation (in 
Method 2 or 4) is on the terms [F^n , F Z £ > Myg only, which are interpo- 

lated from the tabulated input data. The objective of the force implementa- 
tion in Methods 2 and 4 is to force the trim solutions to be accurately corre- 
lated with the wind-tunnel data. Since the trim solution for the rigid case 
depends only on the integrated effects of the pressures, the force correction 
allows the experimental force data to represent global effects of the non- 
linear flow on the rigidly modeled aircraft (i.e., the rigid term in eq . (2)). 

The effect of the pressure implementation in Method 3 is to make correc- 
tions to both the theoretically calculated rigid and elastic terms in the 
right side of equation (2) to account for the local nonlinear flow effects. 

If the experimental force and pressure implementations are applied to the 
FLEXSTAB (Method 4) simultaneously, the effect of the pressure is on the 
elastic increment term only while the effect of the pressure on the rigid term 
is superseded by the force data that are allowed to represent completely the 
global effects or the rigid term in equation (2) . The trim solutions and the 
calculated pressure distribution on a flexible aircraft thus have the combined 
effects from the force and pressure data measured on the rigid wind-tunnel 
model and the results of the theoretically calculated elastic increment. The 
objective of this scheme is to force the trim solutions to be accurately cor- 
related with the wind-tunnel data, and allow the experimental forces and pres- 
sures to represent both the global and the local effects of the nonlinear flow 
while the elastic effects are calculated theoretically. Since the linear 
theory alone cannot treat the nonlinear flow phenomena which are inherent with 
the Orbiter conf iguration , an effort was focused on seeking a correction 
approach in which the linear theory solutions are implemented with experimen- 
tal data. Methods 2 and 3 are investigated before Method 4 in order to find 
an approach by which the linear theory is integrated with less experimental 
data than is Method 4 and yet one capable of predicting the remaining experi- 
mental data correctly. Although Methods 2 and 3 are shown to be inapplicable 
to the Orbiter, it should be noted that their techniques (appendix) are valid 
and may be applied to other aircraft configurations. The application of the 
above correction techniques to a realistic aircraft has never been reported. 
This is probably due to the lack of experimental data and at least partially, 
to the nature of the technique, which requires a time-consuming effort to pre- 
pare the required correction data to make the data compatible with the 
FLEXSTAB system. This report presents the results of computations based on 
the linear theory and/or the linear theory implemented with experimental data. 


7 



Dynamic stability analysis- The linear aerodynamic theory in FLEXSTAB 


(ref. 1) is used to compute the pitch parameters, 
Orbiter in the Mach number range of 0.3 to 1.2. 




and for a rigid 


RESULTS AND DISCUSSION 


Static Stability and Elastic Deformation 

Four analytical methods were used to determine which method gives static 
stability and the lifting pressures comparable with experimental data, and 
thus the one that can be used to estimate correctly the aeroelas ticity effect 
and the aircraft deformation. The results of computations based on the fol- 
lowing methods are presented. All of the results are for the Orbiter config- 
uration depicted in figures 1 and 2. The computation for the rigid Orbiter 
model is compared with the wind-tunnel data for Mach number 0.9 at a = 10° 
and 6=0°. 

Linear theory method- The calculated lifting pressures are shown in 
figure 6. On the fuselage surfaces at x = 13 to 25 m where potential-like 
flow exists due to the straight contour of the fuselage, the calculated pres- 
sure is in good agreement with the experiment. On the nose and tail regions 
of the fuselage, discrepancies of pressures are noted between the calculated 
and the experimental results, but the qualitative trend of the calculated 
pressure distribution is similar to that measured on the wind-tunnel model. 
This indicates that the pressure discrepancies are strongly associated with 
the nonlinear flow occurring on the nose and tail regions of the fuselage. 

The lifting pressures on the wing at various spanwise stations are shown in 
figures 7(a)-(g). It is noted that the agreement between the calculated pres- 
sure and the measurements is fairly good in the midchord region of the 
inboard wing section. The agreement, however, is poor in the region near the 
leading and the trailing edges and also on the outboard wing section. These 
discrepancies are attributable to vortex flow on the leading-edge region and 
on the outboard section of the wing, and to flow separation near the trailing- 
edge region. 

The zero aerodynamic coefficients calculated by using linear theory are 
shown in table 1. The calculated Cl and are comparable to the experi- 

mental results within an order of magnitude, whereas the calculated 
shows a small drag relative to the experimental result. The calculated static 
stability coefficients are shown in table 2. The comparison of the results of 
the linear theory (Method 1) with that of Methods 2 or 4 (which are essen- 
tially interpolated from the wind-tunnel data) shows that the agreement 
between the results from the linear theory and experiments is generally good, 
within an order of magnitude. The agreements in and are 

excellent. Q Q 

Force correction method- This method provides trim solutions for which 
the trim angle is calculated in an iterative manner from the tabulated experi- 
mental data. The lifting pressures are then calculated using the linear 
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theory for the trim angles. The computation results in trim angles a = 9.5° 
and 6 = 7.4°, and the corresponding pressures shown in figures 8 and 9. 

These results are compared with those from the linear theory and experiments 
which specify a = 10° and 6=0°. It shows that the effect of the elevon 
deflection is to increase the pressure on the tail region of the fuselage and 
on the trailing-edge region of the wing. However, this method still does not 
predict accurately the experimental pressures (fig. 9) . It is obvious that 
the linear aerodynamic representation of the actual flow by this method is 
poor and thus Method 2 would not be a good choice for studying the aeroelastic 
effects. 

Pressure correction method - Incorporation of the experimental pressure 
data into FLEXSTAB through the use of equation (1) results in the distribu- 
tions of pressures given in figures 6 and 7. The match of the results of this 
scheme and those of experiment is excellent, indicating that the experimental 
pressures are correctly incorporated into the linear theory solution. The 
aerodynamic derivative results from this method are shown in table 3 and are 
compared with those interpolated from experiments. It is noted that the 
results for Cl o > Clq> Cdq> anc ^ ^Mq are *- n g°°d agreement. However, the dis- 
agreement of c m 0 > C m^> Cl^ is rather contradictory to the fact that these 

parameters are related to the aerodynamic forces that are mathematically inte- 
grated from the experimental pressures and thus should agree with the measured 
forces and moments. Since c m Q ? C m ^, and are associated with the elevon 

it leads to the belief that the disagreement is attributable to the Orbiter 
elevon which is geometrically defined in the computation differently from the 
wind-tunnel model; that is, the trailing-edge region of the wing is smoothed 
and the edge is assumed to be closed in the computation whereas the Orbiter 
actually has a blunt trailing edge with thickness varied from 7.62 to 17.65cm. 
This method is considered inappropriate for use in the aeroelasticity study 
since some of the aerodynamic forces and moments are still not correctly 
predicted . 

Force and pressure correction method - Because of the foregoing results it 
was decided that the linear theories implemented with force and pressure cor- 
rections would be used to compute the aeroelasticity effect and the Orbiter 
structural deformation. Table 4 shows the calculated results of trim solu- 
tions including the aerodynamic coefficients and derivatives for the rigid and 
elastic cases. The comparison of the rigid and elastic cases, both at the 
same flight condition (1 g or 2.5 g) , indicates that the aeroelastic effects 
are negligibly small. However, it is noted that: 

1. A slightly larger elevon trim angle is required in the elastic case 
than in the rigid case to account for the aeroelastic effects. 

2. As the Orbiter is assumed to be flexible, C m is increased but C m 

a o 

is decreased. This results in practically no change between .the rigid and the 
elastic cases at trim condition in the longitudinal stability, C m . 

3. The elevon effectiveness in controlling lift (Cl^) and longitudinal 
stability (C mo ) is decreased when the Orbiter is assumed flexible. 
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4. The aeroelastic effect on lift and drag coefficients is negligibly 
small, although a small amount of loss in lift is evidenced by a small 
decrease in Ct or Ct • 

5. The negligibly small aeroelastic effect on and indicates 

that the dynamic stability study may be confined to the^rigid case only. 

The effects of g loading on the trim conditions as observed from table 4 are 
as follows for both the rigid and elastic cases: 

1. Higher g-load condition requires higher trimming of angles a, 0, 
and 5 . 

2. The zero derivatives Cl q and are not affected by the increase 

of g-loading whereas the derivatives Ct and C m are greatly affected. 

a a 

Thus the Orbiter at high g-load condition loses its effectiveness in lift and 
longitudinal stability control. 

3. The effect of g-load change on C D 0 and C Da is negligibly small. 

4. The derivatives of C^, Cp, and C m with respect to Q and 6 are 
practically independent of g-load condition. 

Figure 10 shows the estimated pressure distributions on the fuselage 
under trim condition in the two flight cases. The pressure increment due to 
aeroelastic effect is negligibly small so that figure 10 actually represents 
the results of both the rigid and elastic models. The results are compared 
with the experimental data for which a and 6 are 10° and 0°, respectively. 
The estimated pressures for the trim 1-g flight case (a = 9.5°, and 6 = 7.6°) 
agrees with the experiment except near the tail (x = 28 m) where the esti- 
mated pressure shows a peak indicating the interference from the elevon 
deflection (6 = 7.6°). On the upstream region of the fuselage ■ (x = 0 to 10 m) 
the pressure varies greatly due to the abrupt change in camberline slopes. 

In the midsection of the fuselage there still exists appreciable pressure 
variation which indicates that the pressure is recovering from the abrupt 
pressure change in the forebody. Figure 10 also shows that the pressure is 
appreciably higher in the 2.5-g flight case than in the 1-g flight case. 

Figure 11 shows the calculated pressure distribution on the wing when the 
Orbiter is in the trim condition. Since the calculated trim angles (a = 9.5° 
and 5 = 7.6°) for the 1-g load condition are not exactly the same as those of 
experiment (a = 10° and 6=0°), the computation shows a pressure rise on 
the elevon area. The pressure rise due to elevon deflection covers only the 
elevon area on the inboard section of the wing. On the outboard section of 
the wing the pressure rise is extended farther upstream of the elevon area. 

The trend of the chordwise pressure distribution on the wing is similar for 
both the 1-g and 2.5-g load conditions. However, in the outboard region of 
the wing where the semispan is larger than 44.20%, the high angle of attack in 
the 2.5-g flight condition gives rise to a sudden pressure increase in the 
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leading-edge region. This pressure rise produces a large pressure differen- 
tial between the leading-edge and midchord region of the wing and thus causes 
total twist on the wing outboard section. On the wing near the 27.1% semispan 
a significant pressure rise occurs at x/ c = 0.25, due to the interference of 
the pressure rise on the adjacent fuselage. 

The results of the estimated fuselage displacements are shown in fig- 
ure 12 for 1-g and 2.5-g load conditions. The fuselage appears to deform in 
the inverted ,l U-shape ,t with maximum displacement occurring at the downstream 
end of the fuselage. It is noticed that a rather sharp inflectional displace- 
ment occurs on the fuselage at x = 19 m and 27 m. The inflection at 
x = 27 m is attributed to the structural discontinuity at that location where 
the cargo bay meets the fuselage. The inflection at the middle region of 
fuselage, x = 19 m, is caused by the first mode bending of the fuselage 
structure. 

Figures 13(a)- (i) shows the estimated wing structural displacements at 
span locations corresponding to those of the pressure distribution shown in 
figure 10. The figures display increasingly larger displacements toward the 
outboard regions of the wing. A very little twist to the wing is observed 
for the 1-g load condition. However, in the 2.5-g maneuver case, a twist 
occurs on the trailing-edge region of the wing, accounting for the pressure 
rise due to the elevon deflection. A twist also occurs toward the leading 
edge of the outboard region of the wing with semispan larger than 61.28%, 
accounting for the sudden pressure rise due to the large angle of attack in 
the 2.5-g flight condition. The twists that occurred in the leading- and 
trailing-edge regions are in the opposite directions and thus result in a 
concave- shaped displacement on the outboard region of the wing. Twist does 
not occur on the leading edge portion of the inboard region of the wing except 
near the extreme inboard section of the wing (27.1% semispan) , where a twist 
occurs near the leading edge due to interference from the pressure rise on the 
adj acent fuselage . 


Dynamic Stability 

The characteristics of the pitch damping of a rigid Orbiter model was 
investigated using the linear theory approach of the FLEXSTAB (ref. 1). 

Figure 14(a) shows the calculated pitch damping represented by and 

in the Mach number range of 0.3 to 1.2. The corresponding + C m ~ are 

shown in figure 14(b) together with the experimental results (ref. 5). The 
following information on the damping characteristics can be drawn from 
figure 14: 

1. The calculated is in good agreement with the experimental 

results . 

2. The negative value of indicates that the Orbiter config- 

uration possesses a stable pitch damping^at subsonic and supersonic speeds. 
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3. As the Mach number is increased from 0.3 to 0.8 the value of 
c n4 + C mQ ^creased, indicating that the damping is increased with 

increasing subsonic Mach number. This increased stable condition is attrib- 
uted to the decrease in both terms C m ? and C^, as shown in figure 14(a). As 

the Mach number is increased to 1.2, C m £ + Cm^ is increased, indicating that 

the Orbiter becomes less stable in pitch damping with increasing supersonic 
speed. This less stable condition is attributed to the increase from 0.85 to 
1.45 of the term 0^. 


CONCLUSIONS 


The FLEXSTAB analysis of the longitudinal stability of a Space Shuttle 
Vehicle-Orbiter model discloses the following. 

1. Experimental aerodynamic force and pressure correction to the linear 
theory is required in the static stability analysis to account for the non- 
linear aerodynamics associated with the Orbiter configuration and thus to give 
an improved prediction of the aeroelastic effect on the stability. 

2. The aeroelastic effect is small. However, it is evident that: 

a. A slightly higher elevon trim angle than that in the rigid case 
is required. 

b. The elevon effectiveness and the longitudinal stability are lower 
for the elastic Orbiter model than for the rigid model. 

c. The aeroelastic effects on the lift coefficients Cl q and 

and on the drag coefficients Cp Q and Cj are negligibly small. 

d. The calculated Orbiter deformation shows an inverted M U-shape M 
fuselage displacement and a concave shape displacement on the 
outboard region of the wing. 

3. The g effect on the trim solutions is evident for both the rigid 
and elastic Orbiter models. The effects are: 

a. Higher g flight condition requires higher trim angles, a, 0, 
and 5 . 

b. Higher g flight condition results in lower a effectiveness on 
the lift control and longitudinal stability. 

c. The drag of the Orbiter and the derivatives (Cp-, Cp-, C^) and 
(Cl^, C D 5 > C m ^) are independent of the g condition. 
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4. The dynamic stability analysis of the pitch damping, C m £ 4- 
rigid model in a Mach number range from 0.3 to 1.2 indicates that: 



for a 


a. The calculated pitch damping is in good agreement with the exper- 
imental results. 

b. The Orbiter possesses a stable pitch damping within the Mach num- 

ber range investigated. However, the pitch damping stability is 
increased with increasing subsonic Mach number, and is decreased 
with increasing supersonic Mach number. This decrease in the 
stability is attributed to the increase toward higher 

positive values with the increasing supersonic Mach number. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, November 7, 1977 
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APPENDIX 


AERODYNAMIC FORCE AND PRESSURE DATA INCORPORATED INTO THE 
LINEAR THEORY SOLUTION - A FLEXSTAB CORRECTION 


A correction method of the FLEXSTAB by incorporating experimental aero- 
dynamic forces and/or lifting pressure data into the solution of the linear 
theory is described. The objective of the method is to amend the linear 
theory solution to account for the nonlinear flow effects. The data to be 
used for such purposes can be those obtained either by experimental or analyt- 
ical methods; however, in the present computation the former approach is used. 
The following describes first the linear theory by which the FLEXSTAB normally 
solves the static aeroelas ticity problem and then introduces the correction 
method to account for the nonlinearities in the calculation of aeroelasticity 
effects. The developments of this correction method described in this appen- 
dix are the summary of those given in reference 1. The essentials of the 
correction method, described in the voluminous FLEXSTAB program are presented, 
and the validity of the techniques are delineated. Some mathematical formu- 
lations in FLEXSTAB are extended to improve the verification of the techniques 
presented in FLEXSTAB. Mathematical terms in equations not essential for 
showing the technique are given only for their physical significance. The 
detailed mathematical expression of these terms can be found in reference 1. 


The Static Trim Problem 


An aircraft in steady reference flight can be described by the following 
classical equations of motion: 


A T 

MU i (Q 1 tan a 1 - R x tan S ^ ) + Mg sin 0 1 = F + F 


XB, 


MUjCRj - P 1 tan a 1 ) - Mg cos e x sin ^ = F^g + F P g 


MU 1 (R 1 tan - Q ^ ) - Mg cos 0-^ cos (j> 1 = F^ + F T 


ZB-, 


(I 


(I zz 1 i yy 1 ^i r i i xz 1 p i q i m xb 1 + m xb 1 
! " I ZZ 1 )P 1 R 1 + I XZ 1 ( ' P 1 2 “ R l 2 ^ _ m yb 1 + m yb 1 


^ t yy, “ I xx 1 - )p i Q i + I xz R i Q i ^zb, + 


M, 


1 


ZB. 


(Al) 


where M denotes the aircraft mass. The components of force and couple, the 
last terms of equation (Al) , produced by propulsion system are assumed to be 
linear functions of the thrust Tj, for example, M^ = T 1 . The present 

study treats the reentry flight problem assuming T 1 = 0. 
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The aerodynamic forces are 
motion a l9 B x , P x , Q 1# R x , and 
and 5^. 




L Tn 


XB “l 

a 


assumed to be linear functions of the aircraft 
the control surface setting, 6e x , 6a l5 



+ 




1 


\ 



4 + 4 p i + 4 R i + 4, Sa i + 4, 5r i 


ZB, 


4 0 + 4/i + 4/> + 4 s /*i 


4, - 4/i + 4 p p i + 4/1 + 4 { / s i + 4 { / r i 
4! ■ 4 0 + 4/i + 4/i + 4,/ e i 


Mi - 


ZB 


4/i + 4/i + 4/i + 4. Sa i + 4. 6r i 

3 P R 6a 6 r 


(A2) 


Equation (Al) contains 12 parameters (Uj , Vj , Wj , Pj, Q 1# R 1 , Gj, T 1 , 

5e, 6a, and 6r) , 6 of which need to be specified; the remaining unspecified 
parameters are solved from the six equations of motion. In the present compu- 
tation, six parameters — U 1# P x , Q x , R x , <i> x , T — are specified in which T = 0 
is assumed for the reentry flight. For a steady pull-up with load factor n, 

Q is calculated by using Q = (n - l)g/U. The equation of motion is solved 
for the trim condition by the FLEXSTAB using Newton’s iteration method. 


The FLEXSTAB has the capability of solving linear and nonlinear trim 
problems : 


Linear case - Equation (Al) is linearized as follows: 

MUj ((*!<»! - RiM + Mge-! = F^ 
MUi ( R i - PjO^) - Mg sin <t> 2 = F^ 
MU 1 (R 1 6i - Qj) - Mg cos 4> 1 = F^ b 


^zz x ~ i yy 1 ^i r i “ I xz 1 p i Q i = 

( - I XX 1 ~ I ZZ 1 )P 1 R 1 + I XZ 1 (P 1 2 ~ R l 2) = 

( ' I YY 1 ~ I XX 1 )P 1 Q 1 + ■’•XZ^l^l = ^Bj j 
The equation (A3) is solved using the linear aerodynamic assumption of (A2) . 


(A3) 
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Nonlinear case - Equation (Al) is solved for two nonlinear cases: 

1. Linear aerodynamics and nonlinear kinematics. In this case the 
linear aerodynamic force representation of (A2) is used in conjunction with 
equation (Al) to formulate the trim problem. 

2. Wind-tunnel aerodynamics and nonlinear kinematics. In this case 
equation (Al) is solved with the aerodynamic force data (nonlinear) obtained 
by experimental or analytical methods. 

The FLEXSTAB yields the following trim solution: the control surface 

setting (Se-p 6a 1 , (Sr-^, the angle of attack [o^ = tan" 1 O^/Uj) ] , the angle of 
sideslip [ftj = tan" 1 (V 2 /V 1 ) ] , and the flight-path, angle (y 2 = Oj - o^). 


The Linear Theory 

Linear aerodynamic equations - Analytical expressions for the aerodynamic 
forces (terms on the right side of equation (Al) are presented. Aerodynamic 
forces for rigid and flexible aircraft are given separately so that the aero- 
elasticity effects can be distinguished. Aeroelastic effects from the propul- 
sion system are not considered. 

The Reference Axis System force and moment vector can be transformed into 
the Body Axis System by a simple sign change for each of the vector components 
since the two systems are colinear but oppositely directed, that is. 


< 4 j i ■ 


(A4) 


where 


and 


^ f b^i ^ f xb’ ^zb’ ^yb’ f yb’ ^xb’ 


[S] 


-1 


-1 


-1 


-1 


-1 


-1 


The aerodynamic forces f*c> are attributed to the local aerodynamic forces at 
each of the n panel segments. They can be expressed as follows: 


{F^) i = 2 [(f>*] T { f A } x 


(A5) 
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where 



■[<>?]“ 


{f A }' 

[<►*] s 

[+J] 

and {f A } x = 

{f A } 




J f n>. 


The matrix [cf)*] is termed the rigid-body mode shape matrix and expressed as 


‘ 1 

0 

0 ’ 

0 

0 

1 

z. 

X 

0 

-X. 

X 

0 

1 

0 

0 

-z. 

X 

Y. 

X 

-Y. 

*- X 

X. 

X 

0 _ 


X ± , Y-^, and are the coordinates (in the Reference Axis System) of the cen- 

ter of pressure at the ith panel or slender body segment relative to the 
center of mass. 

The term {f^}, denotes the vector of the slender, interference, and thin 
body forces expressed in the Reference (Mean) Axis System. Based on the 
linear aerodynamic representation in FLEXSTAB, {f^}, can be written as 

{f A } 1 = q 1 [T fT ][TRANS t ]r' q fc ^^{Cp} 1S ° + f* AM^. 

+ q^T^HT^HT^l^sKCp^ + q x [T fT ] [TRANS^] f' q £ >] 1 {D} 1 (A6) 



where , 

<S } i ■ + {e * } i + { *m j i> + tc P )lnt <A7) 

represents the lifting-pressure coefficients. The steady flow incidence 
can be expressed as: 
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(A8) 


{ Vl " {+«J 4 a, + + {*« r } 5r l + ( *a>“l 

+ { *6>Sl + <*p> | *1 + <* Q > I Q, + %> I K 

Physical significance of the transformation matrices in equation (A6) are 
described as follows: 

[Tf x ] — A matrix that transforms the pressure from the local thin and 
slender-body axis system to the Reference Axis System (eq. 3.5-30, ref. 1). 

[TRANS t ] — A matrix that transforms the thickness induced pressure, that 

is, {Cp}^ S °, to components of forces in the local thin and slender-body axis 
systems (eq. 3.5-31, ref. 1). 

[T T p] — A matrix that transforms forces on the mean camber surface to the 
components expanded in local Body Axis System (eq. 3.5-32, ref. 1). 

[Tpp] — A matrix that transforms lifting pressure q{Cp} on a panel ele- 
ment to the components of force in local thin and slender-body axis 
(eq. 3.5-33, ref. 1). 

The various terms on the right side of equation (A6) can be given the 
following general interpretations: 

1. First term: There are two components; one representing the effect of 

isolated thickness on body forces and the second representing the effect of 
Mach number variations on the forces induced by isolated thickness effects. 

2. Second term: This term accounts for the effects of lifting pressures 

and interference pressure on the body forces. The lifting pressures here 
represent those caused by camber shape, elastic rotation, aircraft motion, and 
control surface deflections. 

3. Third term: This term represents the components of body forces aris- 

ing from leading-edge panel force. 

The second term is of primary interest to the present study because, as 
indicated in equations (A6) to (A8) , it is dependent on the problem motion 
variables. Note that {Cp}, in equation (A7) , is partially computed from the 
operation of the steady aerodynamic influence coefficient matrix [Ap q ] ^ on 
the sum of cambhr incidence vector { ^ 1 » t * ie m °tion incidence vector {^} 1 
involving aircraft rotation, and control surface deflection, and aircraft 
deformation {0*}. 

Linear structural equations- The structural displacements {6} 1 are com- 
puted from the applied aerodynamic force, {Q A }, by 

{6}! = [C] 1 {Q A } 1 (A9) 
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where [C] 1 is termed the unconstrained flexibility matrix. The deformation 
incidence vector {0}* is related to {6} 1 by 

{e*} ! = [p 0 ]{5} 1 (A10) 

A 

The airloads {f-j.^ of aerodynamic load points can be transformed to struc- 
tural node forces by 

{Q A }j = [P t ] T { 4>! (All) 


Substitution from equations (A9) and (All) to equation (A10) yields an expres- 
sion for deformation incidence in terms of aerodynamic segment forces, 


(0 } l ^0T^l^ f T^l 


(A12) 


where 


1 


[p 0 ][c] 1 [p t ] j 


Combined aerodynamic and structural equations - In order to obtain a 
linear expression in the unspecified trim variables, the third term in equa- 
tion (A6) representing the leading-edge correction effects is deleted from the 
aerodynamic representation. Substitution of equation (A7) into equation (A6) 
reduces equation (A6) to the following form: 


(f A } 1 = 


(A13) 


where the term (frp} 1 is defined by those terms in the right side of (A6) 
which post-multiply [Tf T ] after leading-edge correction effects have been 
deleted. Furthermore, (f^i is separated into two parts: 


«$, - [ T TF HFA> i + {f M 


(A14) 


where 


"S -3 1 C'Wi (9 * ) 


(A15) 


and 


fp ?> - q 1 tT TF HT FP H'q 0 'l + {*„),) + {cp 

,S^iso 

/ c l dl 

+ q [TRANS 


V, int 


] 




{ C S} iso + ^ AMn] 


3M 


n 


(A16) 
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Explicit expression of {f^i in terms of {0*} can now be eliminated by sub- 
stituting equations (A12) and (A15) into equation (A14) , that is, 

{f£>l = (M 'J-q 1 [T TF ][A F0 ] 1 [C 0T ] 1 )- 1 {f^}J (A17) 

In order to save the computational time in matrix operations, equation (A17) 
is reduced to the form: 

{fp 1 = (MSI + q 1 [T TF ][D]- 1 [A p0 ] 1 [C 0T ] 1 ){f^}J (A18) 

where the aeroelastic matrix is identified as: 

[Dlj 1 = £" 1 Si ^1 ^ A F6^1 ^0T^1 ^ t tf^ 

The term Ifp ^ in equation (A18) can be identified as 

{fp[ = q 1 [T TF ][A FQ ]({^ c } 1 + {^j) + {f£} lnt + {4> iS ° (A19) 


in which 


<4 } f ° E 5, [TRANS t ] [• q t J ^Cp}j S ° + MM -j j 

E 5l[ T IF 1 f T FP>f < lc' ltC P }int 

The derivation of equation (A18) can be accomplished by using a matrix 
identity, 

[ T Tf^ i 1 f A F0^ 1 + ^ 1 ^ = ^1] ~ ^1 ^ T TF^ A Fe ^ 1 ^0T^ 1^ 

Aerodynamic derivatives for a flexible aircraft - Combining equations (A4) , 
(A5) , (A13) , (A18) , and (A19) leads to a complete expression for the total 
aerodynamic forces in the Body Axis System: 

«£>, - ([° T ] + V G e l [e 0T ] i> <{f X SO + {f T } i nt) + ^V'Vl + 'VP <A20) 

where 


[G t 1 - 2[S]U*][T fT ] 

[G e ] - [6,.] [6J- 1 [T TF ] [A^] 


(A21) 


The longitudinal aerodynamic derivatives can now be obtained by equating equa- 
tions (A2) and (A20) : 
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(A22) 


fid 


XB, 




ZB r 


4 


L rB o. 


= ( [G x ] + q 1 [G 0 ][C 0 T ])({f^}^ S ° + {f^}^ nt ) + q 1 [G Q ]{^ c } 1 


where 


= q^TRANS^M,. ^{Cp}j; SO 

(A23) 

{fM S ° = 5 1 [t tf ][t fp ]M cS 1{c ^} int 





q 1 [G 0 ]{ ^} 


(A24) 





^V { V 


(A25) 



(A26) 


The subscript E refers to an elastic aircraft. 

Aerodynamic derivatives for a rigid aircraft - The flexible aircraft aero- 
dynamic derivatives expressed by equations (A22) to (A26) consist of two 
quantities — the rigid aircraft aerodynamic derivatives and an increment due 
to the static aeroelastici ty . This is accomplished by decomposing the matrix 
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[G e ] into rigid and aeroelastic increment components by using the identity of 
the aeroelastic matrix [DJ” 1 ]: 

[D]^ 1 = M si + qjlT^HDl-ilApgljtC^] (A27) 

For a rigid aircraft the flexibility matrix vanishes; that is, [Cq^] = 0, and 
therefore [D ]“£ 1 = ^ I >]. 

The aeroelastic component of [D]^ 1 can, therefore, be expressed as 
[AD” 1 ] = [D ]*J 1 - fv I >J . It follows that [Gq] of equation (A21) can also be 
separated into two components: one for rigid aircraft component [Gq r ] and the 

other for aeroelastic increment component [Gq E ]; that is, 

G 0 ‘ “W + ^SE 1 (A28) 

where 

^ G 6e^ = [g t 1[ad 1 ] [t tf 1 [a F0 ] x 

For a rigid aircraft the elastic rotation term {0*} vanishes in equa- 
tion (A12) , that is, [Cq T ] = 0, and the aerodynamic forces become 

{F^} 1 = [G T ]([f^ S ° + + q 1 [G 9R ]({^ c } 1 + U M }j) (A29) 

The longitudinal aerodynamic derivatives can be obtained by equating equa- 


q l ^ G 0 R^^c^i 


(A30) continued on next page 


tions (A2) and (A29) : 




XB, 




ZB, 


4 > 


L * B o. 


= [G x ]({f^ S ° + (f x }^ nt ) + 


4 b c 

*4, 


= q l [G 6R ]{ ^ } 


'R 
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XB„ 


F^. 


ZB, 


q l ^ G 6R^Q^ U 


4 


' XB. 


ti 


ZB, 


= q l ^ G 0R ] { ^6 } 


^YB 


•R 


(A30) concluded 


Aeroelastic 'increment to the aerodynamic derivatives - Substitution from 
equation (A28) to equation (A22) separates the aerodynamic derivatives for a 
flexible aircraft into two parts: one for the rigid aircraft and the other 
for the aeroelastic increment; that is. 


(A31) 


1 


> 

1 



^A 


^A 


^A 

'“a 

4* 

— 

f a 

zb a 

4 

+ 

AF^ 

zb a 

4 4b. 

X 

L J 

E 

X 

R 

X 


where the subscript X takes 0, a, Q, or 6 e . The aeroelastic increment to 
the aerodynamic derivatives can be written as 





4i zb 0 


- 5ilS I[5 eT I(( 4>l SO + {f T ,int > + V G 6E 1{ *c } 


E 1 





4F £b„ 



A 

E f 

> 

(A32) continued on i 
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- q l ^ G 0 E-^q} Uj 


"V 

AF^ 

ZB„ 


SM ?i 


qJe' 


r< 


Iaf^ 


ZB . 


q i *- G eE-^6 * 


4m Jb 


Je 1 


(A32) concluded 


Implementation of Experimental Aerodynamic Force Into the Linear Theory 

Force correction- The approach of implementing the experimental aerody- 
namic forces into the linear theory solution of the FLEXSTAB is based on the 
theoretical foundation that aerodynamic derivatives are composed of rigid com- 
ponents and elastic component increments. 


The effect of the experimental aerodynamic force incorporation is on the 
rigid term only in equation (A31) . The elastic increment term on the right 
side of equation (A31) is calculated based on the linear theories represented 
by equation (A32) . 

In this approach, the rigid component of the trim solution is forced to 
represent the wind-tunnel data accurately, and to allow the experimental data 
to represent global effects of the nonlinear flow. The linear theory is to 
use the rigid wind tunnel results as a baseline to estimate the elastic 
effects of the actual airframe configuration. 


The aerodynamic force data are supplied to the FLEXSTAB system in terms 
of C L , C-Q, and C m which are tabulated for various angles of attack for con- 
stant elevon deflections. Since the data for C^, Cp, and C m are measured in 
reference to a Stability Axis System, the data are transformed to a Body Axis 
System in order to be compatible with equation (Al) . The equation of trans- 
formation can be written as follows: 


”4" 


C L 

4 

= [SD S ] 

C D 

4 

W md 

R 

C m 

wm — 


(A33) 


where 


[SD S ] 


is a transformation matrix 


defined in chapter 5 of reference 1. 
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Implementation of Experimental Pressure Into the Linear Theory - 

Pressure Correction 

The effect of the pressure correction is twofold — one on the rigid aero- 
dynamic derivatives and the other on the aeroelastic increment to the aerody- 
namic derivatives. These effects can be demonstrated as follows. The body 
face expressed in the local coordinate system can be written in the form: 

1 ~ ^ T TF^ 1 1 ^ A F0 ]({ *c> + + ^ 1 ^TF^ ^ 1 1 ^TP ^c ^ t C p ^ # 

+ q 1 [T TF ][D]- 1 tA FQ ][C 0T ]{f^}5; SO + {fpi SO (A34) 


In FLEXSTAB, the last term in equation (A34) is omitted so that it can be 
expressed as 

<4 5 1 - IT-rpIfF*}, 


where 

<^>1 " + <* M » + V £ eT Hf T J J S °) + 5 1 [T FP H'q c ')ICp] lnt l 

(A35) 


Equation (A35) transforms the normal force components at the aerodynamic cen- 
troids {F^}-£ expressed in the Reference Axis System to components on the 
local axis system {f^]^. The normal pressure coefficients are defined in 
terms of {F^}^ by 


{AC } = ~ 

P qi 


t T 4p 1{FA) l 


where 



] = 



i 

s zm 2 


S ZM n ’ S ZW X ’ S ZW 2 





Sziy^ an d ^ZW m denote area for nth slender-body segment and mth thin-body 
panel, respectively. Specifically, for a slender body on the plane of sym- 
metry, the area is given as 



ttA. 

i 


where A-^ is the area projection of i segment on x-y plane as shown in 
figure 1. The area of the wing panel, S Z y , is the associate panel area on 
the wing mean surface. m 
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FLEXSTAB allows the user to provide linear corrective inputs to modify 
the pressure computed using the linear theory. Assume that at the ith aero- 
dynamic control point, the plot of ACp b versus angle of attack is as shown 
in figure 15. Suppose for FLEXSTAB, ^ 


while from experiment 


Acf = AC^" + ACjf" a 

1 1 „ J-ra 


An WT . _WT , , _WT 
AC P± = AC P ± + AC Pi ° 

1 Iq lr 


L a 


(A36) 


(A37) 


FLEXSTAB corrects equation (A36) by adding the term 


4c Pl + 4 C p ± “ 
1 0 a 


(A38) 


where 


WT 

AC„ = AC" - AC. 


(A38) 


WT <£ 

4c p - AC - AC^ 


The nonlinear experimental data are, in practice, represented by a linear 
approximation. The data are added on equation (A35) , and the corrected pres- 
sure and aerodynamic force can be expressed as: 


{4C P>c • 17 1 V’A + l4< v + “V 
(f \-VV' 1{4C p> 

- (fA} 1 + 5,[T 4p r l ({AC po } + (AC^lc) 


(A39) 


(A40) 



FLEXSTAB transforms normal pressure to normal forces using the matrix 
] = [T^] - *. The corrected form for equation (A19) becomes 


l4\ - q 1 [T IF HV((' l Vi + { Vl> + {f T }lnt + (f T ,1S ° 


+ q [T TF ][T Fp ]({AC p(i } + {AC^ 


}a) 


(A41) 


and the corrected body force corresponding to (A20) becomes 
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(A42) 


o6 


B J 1 C ^ G T-* + q i *- G 0^^0T-*l^ f T } l c 


The corrected aerodynamic derivatives can be derived by comparing the corre- 
sponding terms in equations (A42) and (A20) . In elastic case. 





([G T ]+q 1 [G 0 ][C 0T ])({f^}^ SO + {f}^ nt ) +q 1 [G 0 ]U c } 1 
+ q x [G t ] [T tf ] [D]" 1 [T Fp ]{AC Po > 


(A43) 


4, 


JE 


q i + q i [ g x H t tf J [6] 1 


[TppHAC } 


(A44) 


In the rigid case: 


4 ( 

4 C 

4. 


4 0 

4 C 

4. 


R 


[G T ]([f^] ±so + {f^} int ) + t G 0R ]{+ c >i 

+ q 1 tG T ][T TF ][T Fp ]{ACp o } 


q lC G 0R ]{,Jj a> + q 1 [G T ][T TF ][T Fp ]{AC Pa } 


(A45) 


(A46) 


Equations (A43) to (A46) gives the corrected aeroelastic increment in the 
following form: 
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- q 1 tG e ][C eT ]({4 J t S ° + <4 lint) + 5l[ G 6E 1W c> 

+ q 1 [G T ][T TF ][AD- 1 ][T Fp ]{ACp a } 



c 


E' 


AF' 


XB r 


4 4. 


e t 


+ ^^FP 


HAC p } 
*a 


(A47) 


(A48) 


The effect of pressure on the aerodynamic derivatives is indicated by the last 
term in each of the above equations, (A43) to (A48) . Note that all other 
derivatives are unaffected. 


Implementation of Experimental Force and Pressure into the Linear Theory - 

Force and Pressure Correction Method 

As the experimental force and pressure are implemented into FLEXSTAB, the 
supplied force data are used only in the rigid part of the computation by 
using equation (A33) ; the pressure data are used only in the computation of 
the aeroelastic increments by using equations (A47) and (A48) . The global and 
local effects are thus achieved by implementing force and pressure data, 
respectively. 
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TABLE 1.- CALCULATED AND EXPERIMENTAL AERODYNAMIC 
COEFFICIENTS FOR A RIGID ORBITER 


Method 

C L 

0 


C 

m 0 

i 

-0.17 

0.001 

0.1 

2 and 4 

-.08 

.08 

.06 

3 

-.09 

.003 

.004 

Experiment 

-.08 

.08 

.06 


(1/deg) 


TABLE 2.- STABILITY COEFFICIENTS FOR A RIGID-ORBITER 
CONFIGURATION (M = 0.9, 1.0 g) 



C L 

a 

C D 

a 

c 

m 

a 

Method 1 

0.04 

0.01 

-0.001 

Methods 2, 4 
and experiment 

.17 

.02 

- . 004 


(1/deg) 



C L 

0 


C 

m s 

Method 1 

0.029 

0.005 

-0.015 

Methods 2, 4, 

.007 

.004 

-.009 

and experiment 


(1/deg) 



C T 

c 

c 


l q 

D o 

m Q 

Method 1 

3.90 

0.73 

-1.99 

Methods 2 and 4 

3.90 

.68 

-1.99 


(1/rad) 



TABLE 3.- AERODYNAMIC DERIVATIVES FROM METHODS 2, 3, AND 4 



C L 

L 0 

C »o 

C 

m o 

Method 3 

-0.09 

0.003 

0.004 

Methods 2, 4, 




and experiment 

-.08 

.08 

.06 


(1/deg) 



C L 

a 

C D 

a 

c 

m 

a 

Method 3 

0.05 

0.02 

-0.002 

Methods 2, 4, 




and experiment 

.17 

.02 

-.004 


(1/deg) 



% 

C °Q 

C 

m Q 

Method 3 

3.89 

0.77 

-0.199 

Methods 2, 4, 




and experiment 

3.90 

.68 

-.199 


(1/rad) 





C 

m 5 

Method 3 

0.029 

0.029 

-0.015 

Methods 2, 4, 




and experiment 

.007 

.004 

-.009 


(1/deg) 
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TABLE 4.- AERODYNAMIC FORCES AND MOMENTS FOR THE RIGID 
AND FLEXIBLE ORBITER AT TRIM CONDITION 



(deg) 


(1/deg) 


(1/rad) 


(1/deg) 



10 r 


EXTENDED 



15 r 


10h 


m 


TYPICAL 

CONTROL 

POINT 



THIN BODY (WING) 


TYPICAL 

CONTROL 

POINT 


TAIL 

32.0575 



WEIGHT: 84163 n 

MOMENTS OF INERTIA (REF. TO BODY AXES) 

l X x = 1.1829 x 108 kg / m 2 
l Y Y = 8.8239 x 10 8 kg/m 2 
l ZZ = 9.0297 x 108 kg/m 2 
l X z = 2.8748x 107 kg / m 2 


— V--i 


J 


10 


15 


20 


25 


30 


35 


x, m 


Figure 1.- Orbiter geometry and structural data. 




Figure 2.- Orbiter wing profile 





(c) Pitch-moment coefficient. 
Figure 4.- Concluded. 

2 i— FUSELAGE AC P vs a (X = 1 1.81 m, = 0.9) 
O W.T. DATA 

EMPIRICAL (LINEARLY FITTED TO 

W.T. DATA 0 < a < 14.5°) 



a 


(a) Fuselage pressure. 


Figure 5.- Experimental lifting pressures representation. 




WING ACp vs a (AT X = 1 1 .653 m, Y = 3.72 m, = 0.9) 
O W.T. DATA 




(b) Wing pressure. 
Figure 5.- Concluded. 


K 

l\ METHOD 1 (LINEAR THEORY) 



0 5 10 15 20 25 30 35 

x, m 


Figure 6.- Lifting pressures on fuselage (Methods 1, 3, 4, and experiment) 
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Figure 8.- I 
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1.0 

AC P 
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0 

- 0.5 

I 

Figure 10.- Est 
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(i) 95.44% semispan. 

Figure 11.- Concluded. 


TRIM CONDITION 



a, deg 

5, deg 

METHOD 4 

ig 

9.5 

7.6 

2.5g 

19.7 

13.7 



\ 


5 10 15 20 25 30 f ' 

TAIL 

32.773 

FUSELAGE X-COORDINATE, m 


Figure 12.- Orbiter fuselage deformation. 









(i) 95.44 % semispan. 


Figure 13.- Concluded. 



MACH NUMBER 


(a) Calculated C m ~ and C^ . 

Figure 14.- Variation of and with Mach number 




(b) Calculated and experimental C mi + . 

a XJ 

Figure 14.- Concluded. 


FLEXSTAB (THEORY) 



a: 


Figure 15.- Variation of AC p 


with angle of attack. 
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